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We presents results of the numerical experiments in favor of the Baez- 
Duarte criterion for the Riemann Hypothesis. We give formulae allow- 
ing calculation of numerical values of the numbers appearing in this 
criterion for arbitrary large k. We present plots of c& for k e (1, 10 9 ). 



1. Introduction. 

In 1997 K. Maslanka pQ proposed a new formula for the zeta Riemann function 
valid on the whole complex plane C except a point s = 1: 
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where coefficients Ak are given by 




This formula was rigorously proved by L. Baez-Duarte in 2003 [2]. In the sub- 
sequent preprint [3] the same author proved the new criterion for the Riemann 
Hypothesis, the journal version of it appeared two years later jl]. The Riemann 
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Hypothesis (RH) states that the nontrivial zeros p of the function ((s) have the 
real part equal 3?(p) = |. Although Riemann did not request it, today it is often 
demanded additionally that zeros on the critical line 3?(s) = \ should be simple. 
Baez-Duarte considered the sequence of numbers Ck defined by: 



'k 



He proved that RH is equivalent to the following rate of decreasing to zero of the 
above sequence: 

Ck = £>(£;- for each e > 0. (2) 
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Furthermore, if e can be put zero, i.e. if Cfc = 0{k~4), then the zeros of £(s) are 
simply. Baez-Duarte also proved in [I] that it is not possible to replace | by | + e. 

Neither in [I] nor in [6] it is explicitly written whether the sequence Ck starts 
from k = or k = 1. However in [1] a few formulas contain k = 0, i.e. summation 
starts from cq. The point is that if we allow k = 0, for which cq = Q/tt 2 , then the 
inversion formula (see e.g. [7j) is fulfilled: 

However I do not see application of the above formula, except the possibility of 
checking some of the statements made in [13] . Furthermore, if the Baez-Duarte 
sequence Ck starts from k = then the following identity holds: 

OO JU oo 



V = p x \^ [ } (a) 
^ k\ ^k\((2k + 2Y y } 



It is an application of the general formal identity: 



J J k\ ^— ' k\ 

k=0 \j=0 vJ/ / fc=0 

where should not increase to fast with k to ensure convergence of series. F] Putting 
here aj = (— l) J 6j gives the usual formula appearing in the finite difference theory 
(see [H] §1): 

k=0 \j=0 VJ/ / ' fe=0 



1 Indeed, collecting on the l.h.s. terms multiplying <Zj we get: Etj (^) ¥ 
a j EfeL/ TET = "./ £n=o 7TT = ".- / ' " and summing over j gives r.h.s. 
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The identity Q can be used to establish the connection with the Riesz criterion for 
RH (original paper [9], discussed in jl], [H]). Riesz has considered the function: 

{x) ~ ^ (k - i)K(2k) -2^k\a2k + 2y 

Unconditionally it can be proved that R(x) = 0(x l / 2+e ), see [H] §14.32. Riesz 
has proved that the Riemann Hypothesis is equivalent to slower increasing of the 
function R(x): 

RH R{x) = O (x 1/4+e ) . (7) 

But from Q we get: 

00 k x 

k=0 

thus the generating function for Ck can be expressed by R(x). In [TD] it is proved, 
that for any real number 5 > —3/2 we have 

R{x) = 0(x &+1 ) ^c k = 0{k 8 ). (9) 



Proof is based on the relation R(k)/k « Ck- 



2. Computer experiments 

The criterion (|2| seemed to be very well suited for the computer verification. 
At the end of [3] Baez-Duarte wrote a sentence "A test for the first up to k = 
1000 shows a very pleasant smooth curve". However for larger values of k the true 
behavior of the sequence turned out to be more complicated: instead of monotonic 
tending to zero there appeared oscillations and Ck changed the sign at first foi0 k = 
19320: C19319 = -1.7870567 x 10~ 13 while C19320 = 9.170232808 . . . x 10" 12 . The next 
sign change is: c 225 26 = 2.2292905301 . . . x 10" 13 but C22527 = -6.5057526 x 10~ 12 . 

To my knowledge the first plot of Ck for k up to 95000 appeared in the book [5] 
published in Polish. The same plot was reproduced in [I]. Data used to make this 
plot consisted of c k calculated every 500-th k — it is very time consuming to get Ck 
directly from ([T]). Indeed, for large j the values of C,[2j + 2) very quickly become 
practically equal to 1, thus the summation of alternating series gives wrong result 
when not performed with sufficient number of digits accuracy. For example, the 
Table I below presents values of the partial sums for C12000 recorded every thousand 
summands (the calculation was performed with precision of 9000 digits). 

Let us remark that the partial sums for n and 12000— n are of the same order. The 
binomial coefficients become very large numbers in the middle and to get accurate 



2 in fact if the sequence starts from k = the first sign change occurs for cq = 6/ir 2 > and 
ci = 6/tt 2 - 90/vr 4 = (6tt 2 - 90)/tt 4 « (-30/tt 4 ) < (more precisely ci = -0.3160113011...) 
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value of Cfc one needs a lot of digits accuracy during the calculation. Maslanka has 
used Mathematica to perform these calculation. Over three years ago I started to 
calculate Ck using the free package PARI/GP p2] developed especially for number 
theoretical purposes and which allows practically arbitrary accuracy arithmetics 
both fixed-point as well as floating-point. I started to calculate consecutive for 
each k with the help of the following script in Pari: 

\p 3500 /* precision set to 3500 digits */ 
allocatemem(250000000) 

range=10000 /* the largest subscript in c_k */ 

denomin=vector (range) ; 

for (n=l, range, denomin [n] =zeta(2*n) ) ; 

default (format, "e22.20") 

{ 

for (k=l, range, c=sum(j=0,k, ( (-1) ~ j ) *binomial (k, j ) /denomin [j+1] ) ; 

write ("c_k.dat" ,k, " ", c)) 

} 



Table I 



n 


2^=0 I L ) V j )((2j+2) 


1000 


8.6575528427959311728 xlO 1492 


2000 


1.0610772171540382076xl0 2346 


3000 


2.6820721693716011525xl0 2928 


4000 


8.4511383022435967124xl0 3314 


5000 


1.8751018390471552047 x 10 3537 


6000 


8.3417729099514988532 x 10 3609 


7000 


1.3393584564622537177 x 10 3537 


8000 


4.2255691511217983562xl0 3314 


9000 


8.9402405645720038417x 10 2927 


10000 


2. 1221544343080764152 xlO 2345 


11000 


7.8705025843599374298 x 10 1491 


12000 


-1.6973092190852083930X10" 7 



The problem I have encountered during these calculations was that it seems to 
be not possible to change accuracy of calculation during running the script (the 
command \p 3500 above). Thus I had to change the precision by hand. It turned 
out that when the precision was to small produced values of were obviously wrong, 
something like ten to the very large power. The rule learned from these examples 
for precision set to make calculations confident was that the number of digits should 
be at least enough to distinguish between 1 and 1 + ^ in the zeta appearing in 
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(lip, i.e. the precision set to calculate Ck should be at least \p = k * log 10 (2). Table 
II presents the real example I have met during calculations: when the precision 
was set to 60000 digits the values of Ck for k between 198000 and 200000 were 
(198000 x log 10 (2) = 59603.93914, 200000 x log 10 (2) = 60205.99913): 



Table II 



k 


Ck 


198000 


-8.1809420017968747912 xl0~ 9 


198500 


-8.1130397250007379108x 10- 9 


199000 


-8.0431163120575296823x 10" 9 


199500 


3.4122583912205353616 xlO 49 


200000 


-1.9276608381598523688x 10 200 



Maslanka kindly send me values of Ck from his calculations up to k = 95000 with 
k jumping in intervals of 500, i.e. k = 500/. Autumn 2005 I have started to continue 
this efforts on the cluster of 8 processors Xeon 2.8 GHz, with 4 GB RAM per node of 
two processors^ with the aim to reach k = 200000 also every 500-th value of k using 
PARI/GP computer algebra system [12]. During last five months of computations 
between 4 and 6 processors I have used to calculate Ck in different intervals of k. 
When these calculations were running I have learned of the paper [4] where explicit 
formulae for Ck in terms of zeros of ((s) were given. Quite recently there appeared 
the paper [6] where the prescription to obtain Ck very quickly were also given. In 
view of these developments there is no need to continue very time consuming cal- 
culations based on the formula 0. The only benefit of these calculation was the 
possibility to compare Ck obtained by means of formulae presented in [I] and [6] 
against those Ck obtained from the generic formula 0. It should be stressed that 
calculations based on ([I]) does not assume the validity of Riemann Hypothesis in 
contrast to formulae presented by Maslanka or below. Using these formulae Ck can 
be calculated very quickly for practically arbitrary k — it is very time consuming 
to calculate Ck without assuming RH. 



3. Explicit formulae. 

The formulae presented in [4] and in [6] expressing Ck directly in terms of the 
zeros of ((s) are essentially the same, they differ in the manner they were derived. 
Maslanka has used the binomial transforms discussed in [H] while Baez-Duarte is 
developing the whole machinery by himself. The formulae of these two authors can 
be written as a sum of two parts: quickly decreasing with k trend Ck and oscillations 

Ck- 



3 because I have used 32-bits version of PARI/GP I was able to use 2 31 bytes =2 GB of RAM 
per process 
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where: 



Cfc 



— y 



Ck = C k + C k 

B(k + l,m) (-l) m (2vr) 2m 



(2tt) 2 ^ Y(2m - 1) C(2m - 1) 



and oscillating part: 



r(fc + i)r(Hg) i ^ Eihll 

Ck -^T(k + l + ^) C'(l ~ P) ~ ~f I 



1 + P 



2 yca-p) 



(10) 



:it 



where it is assumed that zeros of £(s) are simple: C'(p) 7^ an d the sum i s over a ^ 
(i.e. on the positive as well as negative imaginary axis) nontrivial zeros of the C( s )> 
i.e. C(p) — and ^ and 



B(w,z) 



T(w)T(z) 
T(w + z) 



is the Beta function. In fact Baez-Duarte is skipping the trend remarking only that 
it is of the order o(l/k) (Remark 1.6 in |4J). Theoretically the formula for Ck is 
valid in the limit of large k, but surprisingly the numbers produced from the above 



formulae ( 10 ) and ( 11 ) are practically the same as obtained from the generic formula 
([I]) for all k, e.g. already for k = 2 we get c 2 = -0.25699711 from Q, while ([To]) 
and (FLT| give c 2 = —0.256969863 and accuracy increases with k. It suggests that 



the integrals Jk appearing in j3] in the proof of the Theorem 1.5 are decreasing to 
zero rather fast with k. 



First let us consider trend. It can be calculated directly from (10): 



Ck 



1 



;-l) m (27r) 



2m 



(2tt) 2 ^ 2 (k + l)(k + 2)...(k + m)m(m + 1) ... (2m - 2) C(2m - 1) 



(12) 



Table III 



k 


Ck from eq.(12) 


Ck from eq.(13) 


Ck from eq.(|l4) 


1 


-2.60052406393 x 10" 1 


-4.0752814729 x 10 3 


-1.6421193331 x 10 1 


10 


-6.9069591105 x 10" 2 


-2.0455052855 x 10 1 


-1.6421193331 x 10" 1 


10 2 


-1.4804264464 x 10~ 3 


-5.9943727867 x 10~ 3 


-1.6421193331 x 10" 3 


10 3 


-1.6248041420 x 10" 5 


-1.6824923096 x 10" 5 


-1.6421193331 x 10" 5 


10 4 


-1.6403755367 x 10~ 7 


-1.6418022398 x 10~ 7 


-1.6421193331 x 10~ 7 


10 5 


-1.6419448299 x 10" 9 


-1.6420436474 x 10" 9 


-1.6421193331 x 10" 9 


10 6 


-1.6421018816 x 10" 11 


-1.6421113244 x 10" 11 


-1.6421193331 x 10" 11 


10 7 


-1.6421175880 x 10" 13 


-1.6421185279 x 10" 13 


-1.6421193331 x 10" 13 


10 8 


-1.6421191586 x 10" 15 


-1.6421192526 x 10" 15 


-1.6421193331 x 10" 15 


10 9 


-1.6421193157 x 10" 17 


-1.6421193251 x 10" 17 


-1.6421193331 x 10" 17 
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Using this formula I was able to produce every 500-th value of Ck for k = 
500, 1000, . . . 10 9 performing calculations in Pari with 100 digits accuracy in about 



4 hours. For large k I have used following asymptotic expansion of (12): 



1 (W | ^ (13 ) 

k (k + l)(k + 2) V2C(3) 12(fc + 3)C(5) 120(fc + 3)(fc + 4)C(7),T 1 ' 

It can be further simplified to: 

- 1 (2tt) 2 , . 

Ck P2C(3) 1 j 

The comparison of these formulae is given in Table III. 

Now we consider the oscillating part Ck- Since PARI/GP does not have built in 
B(x, y) function, I had to use T(z) functions instead. Because of the fast growth of 



the T(x) function even in PARI/GP it was not possible to pursue with formula (11 ) 
for large k. Namely it crashes for k = 356000 because of overflow. But there is a 
following asymptotic formula (see e.g. [15J, §1.8.7): 



r(x) 

~ X , x — > oo 



T(x + a) 
thus we have 

B(a, x) ~ x~ a r(a) for x large. (15) 

Using it for large k and assuming the Riemann Hypothesis: Pi = \ + ilu pi = 
\ — iji(= 1 — pi) after collecting together in pairs conjugate zeros we get: 



2 



Ck = ~.,3 a i cos Qti lo &( k + !) J - A sin Q7, log(fc + 1) J , (16) 



+ 1=1 
where I have denoted: 




m - V [ I (17) 

A = 



In (16) the decreasing of like k~± is obtained as an overall amplitude of the 
"waves" composed of the cosines and sines with the "frequencies" proportional to 
imaginary parts of the nontrivial zeros of ({s). The coefficients ol\ and j3i decrease 
to zero very fast with I. Namely using the Hadamard product for ((s): 



Uj 2( s -i)r( s /2 + i)llv^ p 



s 
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where C = 0.57721566490153286. . . is the Euler constant, the derivative of ((s) at 
zeros can be computed. Taking into account miraculous simplifications, ppi+ppi = p\ 
and the identity 

r(i + z)r(i-z 



7TZ 



sin (7i z) 



I have obtained that: 



|a ; |oce-^ /4 , |A|oce-^ i/4 . 



(19) 



Because imaginary parts of zeros take large values it suffices to sum in (16) over a 
few first zeros. I have used 10 zeros and the table below gives coefficients a.\ and j3i 



and comparison with (19). 



Table IV 



I 


a t 


A 




1 


2.029173866 xKT 5 


-3.315924256 x 10~ 5 


1.50914 x lO" 5 


2 


-3.333265938 xKT 8 


-1.298336420 x 10~ 7 


6.75315 x 10- 8 


3 


2.886139424 xlO" 9 


-4.153918097 x 10~ 9 


2.94404 x 10" 9 


4 


4.813880001 xHT 11 


-6.332017430 x 10~ n 


4.19039 x 10" 11 


5 


7.546769513 xHT 12 


7.526891498 x 10~ 12 


5.83506 x 10" 12 


6 


6.162524600 xlO" 14 


1.942118979 x 10~ 13 


1.51209 x 10" 13 


7 


-1.578482027 xlO" 14 


1.184829593 x 10" 14 


1.10374 x 10" 14 


8 


-1.07138189 2xl0" 15 


-2.209146437 x 10~ 15 


1.66491 x 10" 15 


9 


9.328038737 xl0~ 19 


-7.472197226 x 10~ 17 


4.22403 x 10" 17 


10 


1.747829093 xl0~ 17 


1.122667624 x 10" 17 


1.05303 x 10" 17 



However already calculations with the first zero 71 = 14.13472514173469 . . . give 
numbers which differ much less than 1% (see / = 1 and / = 2 in the above table) 



from those calculated with larger number of terms in (16) as well as with Ck for 
k < 200000 obtained directly from ([T]) without assuming RH. The plots of for 
k up to 10 9 obtained from these formulae are given in the Fig.l and Fig. 2. In 
Fig.2 there is logarithmic /c-axis and thus the plot has a constant "wavelength", not 
depending on k like on the Fig.l. The envelope is given by: 

2A 

y = __ A = 0.777506276445256 x 10" 5 (20) 

(k + l) d / 4 



and was obtained in the following way: First I have maintained in (16) only the first 
zero p\ = \+ 71: 

c k = — — —3- (aicos f-7ilog(A; + l) ] - ft sin f - Tl log(fc + 1) J J . (21) 



(* + l) 



4 



2 ' ov ' ' \ 2 
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Next I made use of the identity: 
a cos(9) — b sin(9) = A sin((j) — 9) where A = \/ a 2 + b 2 , <p = arctan ^- j (22) 



to obtain: 



5 k = ~ — ~ — rr + P! sin I - i 7l log(fc + 1) ) . ( 23 ) 



(* + l) 



2 



from which (20) follows and numerical value of A is obtained from ax and Pi in 
Table IV. Here0 = — 0.54916(= 56.497°) Let us remark that this value of A agrees 
very well with amplitude reported by Beltraminelli and Merlini |16j . It is interesting 
to note that lhs of the above formula is valid not only for integer k but also for real 
k, thus using the approximation ~ R(k)/k derived in [10J we can write for large 
x: 



R(x) = 2x 1/l (tucos I -7i log(i) I -ftsin I -71 log(i) ] | (2!i 



There is another way of checking accuracy of the above equation (23). Namely 



assuming that (23) is true and denoting by k and k two consecutive zeros of 



Cfc' = 0, c^' = we get for 7l 



2-7T 

71 = io g ((r + 1)/(*' + 1)) (25) 



To make sense, in the latter approach an independent of (23) and relatively fast 
method of calculating is needed. In fact in [1] Baez-Duarte gives among others 
following formula being the transformation of (jT|: 

n=l V 7 

Using this expression I have searched numerically for sign changes of c& up to k = 
10 9 and Table V presents 71 calculated from the above formula (25) together with 



consecutive zeros extrapolated from integer values using the linear approximation 
(let us remind that sin(x) has derivative ±1 at zeros!). The correct value is 71 = 
14.13472514173469379045725 . . .. 

3 

In Fig.2 the plot of k^Ck is presented. The Baez-Duarte criterion requires this 
"wave" to be contained in the strip of parallel lines for all k. The violation of the 

3 

RH would manifest as an increase of the amplitude of the combination k^Ck for 
large k. This point is elaborated in more detail by Maslanka in [B]. Here I will make 
some further comments on this issue. First it should be remarked that the r.h.s. of 



(16) consists of products of three terms: the first depending only on k (the overall 

3 

factor h*), the second depending only on imaginary parts of nontrivial zeros of £ 
(the coefficients ct\ and j3i) and third ingredients depending both on k and I (the 
trigonometric functions). Assume there are some zeros of £ off critical line. We can 



split the sum over zeros p in (11) in two parts: one over zeros on critical line and 
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second over zeros off critical line. This second sum should violate the overall term 
p resen t j n g rs t sum . Let 7 4 -°' ) denote the imaginary parts of the zeros lying 
off critical line ("o" stands for "off"). It is not clear whether asymptotic similar to 



(19) will be valid for zeros off critical line, but it seems to be reasonable to assume 



that it should not differ significantly from (19). Then the contribution to of such 

— _ (°) 

zeros off critical line should contain a factor of the order e 7 * . Because value of 

the imaginary part of the hypothetical zero off critical line should be extremely 

large, perhaps even as large as 10 100 can be expected, the combined contribution to 

Ck coming from the second sum seems to be extremely small, thus to see violation 

of the Baez-Duarte criterion the values of k should be larger than famous Skewes 



number and 



ook something like 10 . Such a big index k should cause that the 



first term in (16) overcame the smallness of the second term depending only on 7^ 

Table V 



(0) 





extrapolated zero k 


71 from eq. ( 25 ) 


1 
1 


lyoiy.uiynoi 


U.OOoy4oOzlOOolzoO ( 


Z 


ZZOZO.UoololZ 


4U.y loUlOl 1 41o04oUoU 


Q 
O 


A 1 Sfifi 0707/1 1 fi 
41o0o.y ( U 1 41o 


lU.loOO ( OloU ( oUo4oll 


A 

4 


ouuy4. 00 11000 


1 1 .00 / 44oozy { yuuozou 





yoo ( o.zoi4ouy 


1 9 7/i 7/i /IT 1 fi^n^oo^nn^ 
1Z. ( 4 1 441 loouoyzouuo 


6 


149320.1629630 


15.05785832239914317 


7 


236817.0977574 


13.62376898016027559 


8 


366000.4553802 


14.43265811387491946 


9 


573460.6753253 


13.99205468779138406 


10 


891841.3774543 


14.22828231751450665 


11 


1393469.2943691 


14.07955314904129825 


12 


2173554.0482344 


14.13327453400666383 


13 


3387835.5708183 


14.15682029348025265 


14 


5283842.8916393 


14.13659885512530467 


15 


8247263.1465316 


14.11229699162269810 


16 


12864372.4128156 


14.13284964864925287 


17 


20052822.4883780 


14.15424390183719722 


18 


31271608.9210745 


14.14046804681939058 


19 


48805962.1935733 


14.11501826908551253 


20 


76145459.4891092 


14.12608973663948854 


21 


118689214.0783221 


14.15568694594979885 


22 


185076299.0519995 


14.14304477168657435 


23 


288851950.0033488 


14.11488330063446422 


24 


450663149.7663923 


14.12567920473796440 


25 


702517003.8056590 


14.15292988626078933 



The plot in Fig. 2 is a perfect sine of one wavelength thus it gives visual justifi- 
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cation of the above statement that c is determined in fact by the first zero 71. The 
same phenomenon was mentioned by Maslanka in [6]. 



4. Final remarks 



The formula (16) together with a few coefficients ol\ and Pi taken from Table IV 
allows to compute values of Ck for arbitrary large k. Other criteria for RH, like the 
value of the de Bruijn-Newman constant [17], are vulnerable to the Lehmer pairs 
of zeros of ({s). It is hard to see the reason for violation of the inequality \c^\ < 

3 

const /c~4. I have checked, that at the first Lehmer pair P6709 = 0.5 + 7005.062866172 
and p 67 io = 0.5 + 7005.1005646z the derivative has value £'(0.5 + 7005.06286617i) = 
3.2229849698 + 0.74179951875? and similar value for second zero, thus there is no 
chance to get values of violating ^ in this way. It seems to be an open problem 
how to connect the value of the largest k for which \ck\ < const/c" to the number of 
zeros lying on the critical line. Let us mention that for the Li's criterion [18] which 
states that if the numbers 

fulfill A n > for each n then RH is true it is known that if the first n Li's constants 
X n are positive then every zero p of ((s) with \^sp\ < y/n lies on the critical line 
3fy = § P5]. 

After a few months of computer experiments with Ck I believe Baez-Duarte se- 
quence is one the most important and mysterious sequences in the whole mathe- 
matics. 
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Fig. 2: The plot of k 3 / 4 c k for k G (10 6 , 10 9 ). 



